Draft version June 11, 2009 

Preprint typeset using I^T^X style cmulateapj v. 08/13/06 



O 

o 

(N 
3 



O 

u 

43 

9 L1, 
6 

& 

(N 
> 
oo 

On 

t> 

(N 
cn 

o 

&\ 
O 



13 



THE DUSTY NUCLEAR TORUS IN NGC4151: CONSTRAINTS FROM GEMINI NEAR-INFRARED 
INTEGRAL FIELD SPECTROGRAPH OBSERVATIONS 

ROGEMAR A. RlFFEL, THAISA STORCHI-BERGMANN 
Universidade Federal do Rio Grande do Sul, IF, CP 15051, Porto Alegre 91501-970, RS, Brazil 

AND 

Peter J. McGregor 

Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia 

Draft version June 11, 2009 

ABSTRACT 

We have used a near-infrared nuclear spectrum (covering the Z, J, H and K bands) of the nucleus of 
NGC4151 obtained with the Gemini Near-infrared Integral Field Spectrograph (NIFS) and adaptive 
optics, to isolate and constrain the properties of a near-IR unresolved nuclear source whose spectral 
signature is clearly present in our data. The near-IR spectrum was combined with an optical spectrum 
obtained with the Space Telescope Imaging Spectrograph which was used to constrain the contribution 
of a power-law component. After subtraction of the power-law component, the near-IR continuum is 
well fitted by a blackbody function, with T = 1285 ± 50 K, which dominates the nuclear spectrum 
- within an aperture of radius 0'.'3 - in the near-IR. We attribute the blackbody component to 
emission by a dusty structure, with hot dust mass Mhd = (6.9 ± 1.5) x 10~ 4 M Q , not resolved 
by our observations, which provide only an upper limit for its distance from the nucleus of 4pc. 
If the reddening derived for the narrow-line region also applies to the near-IR source, we obtain a 
temperature T = 1360 ± 50 K and a mass M H d = (3.1 ± 0.7) x 1O~ 4 M for the hot dust. This 
structure may be the inner wall of the dusty torus postulated by the Unified Model or the inner part 
of a dusty wind originating in the accretion disk. 

Subject headings: galaxies: individual (NGC 4151) - galaxies: active - galaxies: nuclei - galaxies: ISM 



1. INTRODUCTION 

The presence of a very compact infrared source in 
the nucleus of the Seyfert galax y NGC 4151 was firs t 
suggested by the observations of iPenston et al.l (|1974f h 
who reported photometric variations of the galaxy nu- 
cleus both in the optical and in the near-infrared (here- 
after nea r-IR). A close inspe ction of the light curves 
shown by IPenston et al.l (|1974l ) suggests a delay between 
the near-IR and o p tical variations, recently confirmed 
by iMinezaki et al.l (|2004l ). who could better quantify 
its value as 48lg days. These authors concluded that 
the K— band emission is dominated by thermal radia- 
tion from hot dust located at «0.04pc from the nu- 
cle us. A similar interpretat ion was previously proposed 
bv lRieke fc Lebofskvl (p8lT i. via models fitted to a large 
aperture (~8") nuclear optical and near-IR spectrum, 
concluding that in the near-IR the emission is dominated 
by thermal reradiation by hot dust with temperature 
Thd = 1300- 1500 K 

Mun dell et al.l (|2003l ) have reported neutral hydrogen 
absorption in VLA radio observations with subparsec 
linear resolution, close to the systemic velocity of the 
galaxy, interpreted as due to a clumpy gas layer suggested 
to be assoc iated to a circumnucl ear obscuring torus. Ac- 
cording to iMundell et al.l ((2003) , this absorption is ob- 
served along the line of sight towards the first component 
of the radio counterjet, at <0.1pc from the nucleus, thus 
consiste nt with the estimated in ner radius of ~0.04pc de- 
rived by IMinezaki et al. (2004]). The interpretation put 
forth bv IMundell et al.l ( 2003h - that the obscuring torus 



is not in front of the active galactic nucleus (AGN) - 
is necessary in order to reconcile the radio observations 
with the UV/optical observations which show a blue nu- 
clear continuum and broad emission lines, indicating that 
both the nuclear source and the broad-line region are vis- 
ible. 

The presence of a dust y torus around the nucleus of 
NGC 4151 was disputed bv lSwain et~afl (|2003D on the ba- 
sis of interferometric observations at 2.2/^m using the two 
10 m Keck telescopes. These authors found a marginally 
resolved nuclear source <0.1 pc in diameter - consi stent 
with the observations of both IMundell et al.l (|2003l ) and 
IMinezaki et al.l (|2004f ). but favor an origin for this emis- 
sion in thermal gas from the accretion disk, arguing that 
previous nuclear spectra of NGC 4151 did not show an 
"infrared bump " , in apparent cont r adicti on with the pre- 
vious study bv iRieke fe Lebofskvl (|198lD . Such a bump 
is expected if a dusty torus is present, due to reprocess- 
ing of UV/optical radiat i on from the central engine (e.g. 
Rieke fc Lebofskvl 119811 : iBarvainid 119871 : ISanders et all 



1989). 
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The alt e rnativ e interpretation proposed by 
iSwain et al.l (2003) seems nevertheless to be in line 
with recent results presented by iKraemer et al.l (|2008h , 
who studied the distribution of [Oiii]A5007 and 
[On]A3727 emission in the Narrow Line Region (NLR) 
of NGC 4151. These authors pointed out that there does 
not seem to be a region of complete shadow around the 
apex of the bicone, as significant ionized gas emission is 
found there, although with a lower ioniza tion parameter, 
consis tent with a weaker ionizing flux. IKraemer et al.l 
(2008) propose that the attenuation of the nuclear radi- 
ation could be due to low-ionization absorbers swept out 



2 



Riffel, Storchi-Bergmann & McGregor 



by a wind originating in the accretion disk, suggesting 
that there is no need for a "classical, doughnut-shaped" 
torus in order to explain the observed light distributions, 
in agreement with rece nt modelling for the obscuring 
region around an AGN (|Elitzurll200a ). 

In this work, we present new observations of the nu- 
cleus of NGC4151 with the Gemini Near infrared Inte- 
gral Field Spectrograph (NIFS) at an angular resolution 
of ft' 12 - corresponding to a spatial resolution of a few 
pc at the galaxy - which is much improved relative to 
previous spectroscopic studies. Our observations clearly 
reveal the presence of an unresolved nuclear source in 
the near-IR. We h ave already used the NIFS da ta to 
map the excitation ([Storchi-Bergmann et al.ll2009L here- 
after Paper I) and kinematics (|Sim5es Lopes et al.ll2~009f) 
of the NLR. In the present paper, we use the nuclear 
spectrum, combined with surrounding extranuclear spec- 
tra and a Space Telescope Imaging Spectrograph (STIS) 
optical spectrum, in order to isolate the emission of the 
infrared nuclear source, which is consistent with that pre- 
dicted for the torus pos tulated in the Unified Model of 
AGNs (|Antonuccilll993l) . 

We adopt a distance to NGC4151 of 13.3 Mpc, cor- 
responding to a sca le at the galaxy of 65 pc arcsec -1 
(jMundell et al.ll2003l ). This paper is organized as follows. 
In section 2, we present a summary of the observations, 
in section 3 we present the results and in section 4 we 
present and discuss our conclusions. 

2. OBSERVATIONS 

Two-dimensional spectroscopic data were obtained on 
the Gemini North tele scope with the instrument NIFS 
(|McGregor et aLlfeOOSf) operating with the adaptive op- 
tics module ALTAIR on the nights of December 12 and 
13, 2006. The observations covered the spectral bands 
Z, J, H and K, at spectral resolutions of 4990, 6040, 
5290 and 5290, respectively, resulting in a wavelength 
coverage from 0.94/im to 2.40/im. More details of the 
observations and reduct i on pr ocedures can be found in 
IStorchi-Bergmann et all (|2009[) . 

In the present paper we concentrate on the analy- 
sis of the nuclear continuum. The spatial profiles of 
the nuclear region of the galaxy in the K and J bands 
are shown in Fig[TJ in comparison with stellar profiles. 
The nuclear source - which has a full width at half 
maximuum (FWHM) of 0'.'12±0'.'02 in the K-band and 
0'.'16±0'.'02 in the J-band - corresponding to 8.0±1.3pc 
and 10.4±1.3pc at the galaxy, respectively - is unre- 
solved by our observations, in spite of the fact that their 
spatial profiles are marginally b roader than those of the 
stars in Fig.[T] As discussed in [Storchi-Berg mann et al.l 
(2009), this is due to the much shorter exposures of the 
stellar observations. 

In order to include as much as possible of the nuclear 
source, we have chosen an extraction aperture of ft' 3 ra- 
dius which includes the nuclear flux down to 10 % of the 
peak value in the K-band, as illustrated by the dashed 
lines in Fig.[TJ 

3. RESULTS 

The nuclear spectrum is shown in the Fig[2l together 
with extranuclear spectra obtained through similar 0'.'3 
apertures at distances of ft' 9 NE, ft/9 SW, 1'.'8 NE, and 
l'.'8SW from the nucleus, normalized to the flux of the 
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Fig. 1. — Spatial profiles for the galaxy and of a star at J and K 
bands. The vertical dashed lines mark the extraction aperture of 
the nuclear spectrum (0'.'3 radius). 




Fig. 2. — Nuclear spectrum (top) obtained for a circular aperture 
with 0'.'3 radius and extra-nuclear spectra at the positions 0"9NE, 
tf'9SW, L"8NE and l'.'SSW, normalized to the flux at the nucleus 
at A fa l^tm.The spectra are shown in arbitrary flux units and are 
shifted by a constant for clarity. 

nuclear spectrum at A «s 1/im. Fluxes of more than 50 
emission lines have been measured and listed in Table 1 
of Paper I for the nuclear and two extranuclear spectra. 

From Fig. [21 as well as from spectra shown in Pa- 
per I, we observe that the extranuclear spectra are 
blue at all locations around the nucleus. As this blue 
continuum is observed up to at least 2 arcsec from 
the nucleus, it cannot be due to the AGN and can 
thus be attributed to a young stellar population sur- 
rounding the nucleus. Absorption features from stel- 
lar CO bands are indeed observed in the K band, 
marked by vertical dotted lines in Fig.O Only right 
at the nucleus the continuum is clearly very red. Sim- 
ilar red nuclear continua have been observed for other 
AGNs and attributed to emission from hot dust, pos- 
sibly from a dusty torus around the supermassive 
black hole, with tem perature close to the dust sub- 
limat i on temperature dBarvainisI Il987t iMarco fc Alloinl 
fl998l [20001: iRodriguez Ardila C ontini & Viegas 2005bT 
iRodrfguez-Ardila k. Mazzalavll2006l : iRiffel et al.ll2009h . 

3.1. Nuclear optical continuum 

In order to better constrain the nuclear spectral energy 
distribution (SED), we have used an optical STIS spec- 
trum of the nucleus extracted within a similar aperture 
to that of our near-IR spectrum, obtained from the atlas 
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Fig. 3. — The optical+near-IR nuclear spectrum of NGC4151, together with a fit (continous line) considering the contribution of a 
power-law (dotted line) plus a black body function (dashed line). 



of AGN spectra compiled bv lSpinelli et all dl)06). This 
spectrum was extracted within an aperture of 0'.'2x0'.'l 
and covers the spectral range 0.3-1 /im. We have used 
the continuum region from 9780 to 9820 A to normal- 
ize the optical spectrum to the same flux as our near-IR 
spectrum (which extends down to the Z-band) as there 
is a flux difference of « 30% between the STIS spectrum 
and our nuclear spectrum. The small difference between 
the apertures is not important because the stellar contri- 
bution, when compared to the nuclear one, is small, as 
discussed below. We show the optical plus near-IR spec- 
trum in Fig.[3l The optical continuum is very blue and 
can be well reproduced by a power-law up to w 1 //m, 
where there is a change in slope and a red component 
begins to appear in the near-IR. Although the near-IR 
emission is dominated by this component, the contribu- 
tion of the power-law in the near-IR is not negligible in 
the Z, J and H bands, as discussed below. 

3.2. Modelling the spectral energy distribution 

We have modelled the optical+near-IR nuclear con- 
tinuum of NGC4151 by the sum of two components: 
a power law, which seems to dominate in the optical, 
plus a blackbody function, which dominates in the near- 
IR, as illustrated in Fig. [3] The fit was performed using 
the nfitld task from the STSDAS iraf package, which 
allows the selection of limited spectral intervals for the 
fit. We have excluded the spectral interval between 0.6 
and 0.9/im from the fit as we have estimated that the 
contribution of the Paschen continuum in this interval 
(« 1.3 x 10~ 15 ergscm- 2 s- 1 A- 1 at 8200A), is not neg- 
ligible when compared to the total flux observed in the 
continuum. Moreover, the quality of the optical spec- 
trum in its red end is quite poor, and should not be 



used to constrain the SED. Excluding the above spec- 
tral region, the best fit was obtained for a power law 
F\ oc A -2 - 25 and a blackbody function of temperature 
T = 1285 ± 50 K, which peaks at « 2.3 ^m. The fit of 
the two functions combined is shown in Fig. [3] as a con- 
tinuous line, while the black body component is shown 
as a dashed line and the power law as dotted line. 

We have next evaluated the effect of the underlying 
stellar population in the fit above. A stellar popula- 
tion spectrum was extracted from our near-IR datacubes 
within a ring around the nucleus with inner radius of 0'.'6 
and outer radius of 0'.'9. We have then normalized the 
flux of this extranuclear spectrum to 15% of the peak nu- 
clear flux in the K band, which is the maximum contri- 
bution of the stellar population estimated from a spatial 
profile through the nucleus. Under the assumption that 
the underlying stellar population in the nuclear spectrum 
is the same as the one from the ring (which is supported 
by other spectra from the nuclear vicinity) , we have then 
subtracted the extranuclear spectrum from the nuclear 
one. As we do not have an extranuclear optical spec- 
trum, we repeated the fit for the near-IR only, comparing 
the temperature of the uncorrected to the corrected spec- 
trum, concluding that, within the uncertainties («50K), 
the blackbody temperature is the same. 

The reddening may also affect the results of the fit. 
An estimate of the reddening in the near-IR can be ob- 
tained from the emission lines of the narrow-line region. 
We have used the [Fen]A1.2570^m/A1.6440^m line ra- 
tio to calculate the reddening (as described in Paper I, 
but using the present nuclear spectrum, for an aperture 
of radius 0'/3), resulting in E(B-V)~ 0.6. It is not clear 
that this reddening applies to the optical continuum. If it 
does, it would result in a very steep power law for the cor- 
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rected nuclear continuum : using the extinction curve of 
iRieke fc Lebofskvl (|1985h . we obtain Fx oc A~ 4 5 . Never- 
theless, the presence of broad lines in the spectra and the 
already very blue continuum suggests that the reddening 
of the optical spectrum is much smaller; we may be seing 
the nuclear source through a low extinction window to 
the nucleus. This may apply also to the near-IR nuclear 
source, but, even if the extinction applies, the correc- 
tion in the near-IR is small. A reddening correction of 
only the near-IR spectrum for E(B-V) = 0.6, keeping the 
power-law index the same, results in a revised blackbody 
temperature of T = 1360 ± 50 K. 

3.3. The hot dust mass 

Following [BarvainisJ (jl987l ) and iRiffel eFax] (|2009h we 
can estimate the mass of the hot dust producing the 
unresolved near infrared emission. Firstly, we obtain 
the infrared spectral luminosity of each dust grain, in 
ergs^ 1 Hz~ 1 , assuming that the dust is composed of 
graphite grains, by 

(1) 



L 



gr _ 



4ir 2 a 2 Q„B„{T gI ), 



where a is the grain radius, Q v is the absorption efficiency 
and B u (T gI ) is its spectral distribution assumed to be a 
Planck function with temperature T gr . 

The number of dust grains (A^hd) can be estimated 
from the ratio between the total luminosity of the hot 
dust (L^ B ) and the total luminosity of one dust grain 

(X? 1 ) as ./Vhd = was obtained by integrating 

the flux under the Planck function fitted to the nuclear 
spectrum and adopting a distance to NGC4151 of d = 
13.3 Mpc, while L? r was obtained by the intregration of 
equation[T] for a temperature of T gr = 1285 K. Finally, 
the mass of the emitting hot dust is obtained by: 



4-7T n 

Mn D w —a AWPgr- 



(2) 



Assuming a graphite d ensity of p gr = 2.26 gem 3 
(jGranato fe Danesel Il994h . we obtain M H d = (6.9 ± 
1.5) x 10- 4 M Q . 

If we assume that the near-IR nuclear source is subject 
to reddeninng, we should then use the temperature T gr = 
1360 K, which will lead to a smaller dust mass of Mhd = 
(3.1 ± 0.7) x 10- 4 M Q . 

In table[IJ we present a comparison of the mass of hot 
dust obtained here for NGC4151 with the ones derived 
for other AGNs in previous works. This comparison 
reveals that the hot dust mass for NG C4151 is similar 
to the ones obtained fo r NGC 1068 bv iMarco fc Alloinl 
pOOOh . for N GC 4593 bylSantos-Lleo et all (119951) and for 
NGC 1566 bv lBaribaud et all (|1992l ) and is smaller than 
the masses derived for other galaxies such as NGC 7582, 
Mr k l239, Mrk 76 6, NG C 7469, NGC 3783 and Fairall9 
bv IRiffel et al l (I2009D. iRodriguez-Ardila fc Mazzalavl 
( 2006h. iRodriguez-Ardila. Contini fc Vie gad 

(l2005bf)lMarco fc Alloinl <l!998Tl iGlassI (11992ft and 
IClavel. Wamsteker fc Glass! (|1989l ). respectively. 



4. DISCUSSION AND CONCLUSIONS 

As discussed in the Introduction, ISwain et al.l (|2003l ) 
observed NGC 4151 at 2.2 /im and have marginally re- 
solved a nuclear source <0.1pc in diameter, conclud- 
ing that the emission was due to thermal gas from 



the central accretion disk instead of a dusty torus, as 
they claimed that previous observations did not show 
the bump. Our near-IR nuclear spectrum (Fig. ^ - 
clearly shows this bump and favours emission by hot 
dust as the dominant mechanism for its origin, in agree- 



(Rieke & Lebofskvl 


1981). A more recent near-IR spec- 


trum obtained by 
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also suggests the presence of the bump. 

The modern view of the torus {Elitzur 2008) is not of a 
static doughnut shape structure but of a clumpy wind of 
dusty and optically thick clouds originating in the black- 
hole accretion disk, a view which has re cently been fa- 
vored also in the case of NGC 4151 by iKraemer et al.l 
(|2008ft . a nd which is in partia l agreement with the propo- 
sition bv ISwain et alj (|2003l ). Being clumpy, the dusty 
wind allows the view of the nuclear source through low- 
extinction holes. The presence of such holes is also sup- 
ported by the observation of ionized gas close to the apex 
of the bicone, beyond its presumed walls. This more 
updated view is also consistent with our data and re- 
sults from previous papers discussed above, if we con- 
sider that the near-IR nuclear source can be identified 
with the hottest part of the dusty wind, which is the 
one closest to the nucleus and which emits most of the 
near-IR radiation. According to our data, the temper- 
ature of this dust is T = 1285 ± 50 K and its mass is 
M HD = (6.9 ± 1.5) x 10 _4 M©, obtained by fitting the 
optical and near-IR spectrum. The observation of the 
nuclear continuum and of the broad emission lines sug- 
gest that we are looking to the nuclear source through a 
low extinction window, in spite of the fact that the esti- 
mated reddening from the emission lines of the narrow- 
line region is E(B-V) =0.6. If we assume that this red- 
dening is also affecting the nuclear source, we obtain 
T = 1360 ± 50 K for the dust temperature and mass 
M HD = (3.1 ± 0.7) x 10- 4 M o . The temperatures ob- 
tained are in good agreement with those previously ob- 
tained by IRieke fc Lebofskvl |l981). Comparison be- 
tween our near-IR flux values with theirs shows that the 
nucleus was 3-4 times brighter than in the epoch of our 
observations, suggesting that the temperature of the hot 
dust does not depend on the AGN luminosity. 

Regarding the location of the hot dust, our spatial 
resolution allows us to put only an upper limit to the 
radius of the infrared source of 4pc (half the FWHM 
of the nuclear source). Sim i lar sc ales have been probed 
bv lStorchi-Bergmann et all (|2005[ ). who found a heavily 
obscured starburst closer than 9pc from the nucleus in 
NGC 1097 and arg ued that it could b e located in the out- 
skirts of the torus. iJaffe et all (|2004h using mid-IR inter- 
feromentric observations constrained the dust emission in 
NGC 1068 to originate in a region of only 1 pc diameter. 
In the case of N GC 4151, optica l and n ear-IR monitoring 
observations of iMinezaki et all (|2004l ) provide a better 
constraint on the location of the inner radius of the torus 
at «0.04pc from the nucleus. 

The NLR of NGC 4151 has an approximate biconical 
morphology with projected opening angle of ~75° and 
with its axis oriented along a PA -60/240° (NE-SW). 
The SW side of the bicone is tilted toward us, making an 
angle of 45° with the line-of-sight, according to Da s et all 
(2005), who gives also an opening angle of the cone of 



The torus in NGC4151. 
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TABLE 1 
Masses of hot dust found in AGNs 



Galaxy 



Mud (Mq) Reference 



NGC4151 
NGC 7582 
Mrk 1239 
Mrk 766 
NGC 1068 
NGC 7469 
NGC 4593 
NGC 3783 
NGC 1566 
Fairall 9 



6.9 x 
2.8 x 
2.7 x 
2.1 x 

1.1 x 

5.2 x 
5.0 x 
2.5 x 
7.0 x 
2.0 x 



10-"^" 
10~ 3 

10- 2 

10~ 3 
10~ 3 
10~ 2 



10~ 4 
10~ 3 
10" 4 
10" 2 



This work (no extinction, see text) 
Riffel et al. (2009) 

Rodriguez- Ardila fe Mazzalav (2006) 
Rodrigucz-Ardila J _Corrtini_& i V - iegas (2005b) 
Marco fc Alloin ( 20001 
Marco fc Alloin (19981 
Santos-Lleo ct al. (1995) 
Glass (19921 
Baribaud et al. (19921 
Clavel, Wamstcker fc Glass 



66°, which puts our line-of-sight just outside the cone 
wall. If the torus is oriented perpendicularly to the cone, 
we are looking down into the torus hole. This allows 
us to constrain the half-height of the torus to less than 
0.04 pc for a "doughnut shape" opaque torus. This torus 
should extend to at least 0.1 pc from the nucleus, in or- 
der to contain the n eutral hydrogen clouds observed by 
iMundell et al.1 (|2003f ). As proposed by these authors, the 
NE part of the radio jet would be behind the torus, while 
the SW part would be in front of it. 

We thank the referee Dr. George Rieke for valuable 
suggestions which helped to improve the present pa- 
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